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Abstract. We present mid-infrared (2-12pm) spectra of the microquasar SS433 obtained with the Infrared Space Observatory 
(spectroscopic mode of ISOPHOT and ISOCAM). We compare them to the spectra of four Wolf-Rayet stars: WR78, WR 134, 
WR 136, and WR 147 in the same wavelength range. The mid-infrared spectrum of SS 433 mainly shows H i and He i emission 
lines and is very similar to the spectrum of WR 147, a WN8(h)+B0.5V binary. The 2-12yum continuum emission of SS433 
corresponds to optically thin and partially optically thick free-free emission, from which we calculate a mass loss rate of 
2 - 3 x 10~ 4 M q yr -1 if the wind is homogeneous and a third of these values if it is clumped. This is consistent with a strong 
stellar wind from a WN star. However, following recent studies concluding that the mass donor star of SS433 is not a Wolf- 
Rayet star, we propose that this strong wind out flows from a geometrically thick envelope of material that surrounds the 
compact object like a stellar atmosphere, imitating the Wolf-Rayet phenomenon. This wind could also wrap the mass donor 
star, and at larger distances (~ 40 AU), it might form a dust envelope from which the thermal emission, detected with ISOPHOT 
at 25 pm and 60 pm, would originate. This wind also probably feeds the material that is ejected in the orbital plane of the binary 
system and that forms the equatorial outflow detected in radio at distances > 100 AU. 

Key words. (Stars:) circumstellar matter - Stars: individual: SS 433 - Stars: Wolf-Rayet - Stars: winds, outflows - Infrared: 
stars - X-rays: binaries 


1. Introduction 


SS 433 was the first microquasar (an X-ray binary with rela¬ 
tivistic jets) discovered in the 1970’s, when that designation 
had not even been invented yet. It was first discovered as a 
star with strong Her emission lines and so was included in the 
Stephenson and Sanduleak catalog ( IStenhenson & Sanduleakl 
1977 ) as object number 433Ht was later associated to a variable 
X-ray and radio source dZ witter et alj fl989). Detailed study of 
its visible spectrum revealed that this source is unique._For 
a com p lete review of SS433, we refer t he rea der to IMareon 
dl984l) .iZwitter et al. Il989i). Vermeulen! Il996l) . and to the in- 
troduction of lGies et aid ( 2002hl) . 

There are two sets of optical emission lines in the spec¬ 
trum of SS433: the first set corresponds to the so-called “sta¬ 
tionary” lines dMargon et alJ 1 9791) showing normal Doppler 
shift movements, including the strong Ha lines, with a pe¬ 
riod of 13 days. The second set regroups the “moving” lines, 
those that are less intense (by about 1/3) but that show huge 
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Doppler shifts (+50 000 Km.s 1 and -35 000 Km.s 1 at the 
maximum elongation) corresponding to relativistic velo cities 
with a period of about 162 days dMargon et ali 11979 ). The 
“kinematic model” is generally the one chosen to explain these 
unusual observations: SS 433 is an X-ray binary system orbit¬ 
ing in 13 days and emitting relativistic jets where the mov¬ 
ing lines are formed. These jets undergo a precession move¬ 
ment in 162 days. The parameters of thi s sy stem were re¬ 
cently re-calculated by lEikenberrv et al ■ d200lb . who found: 
F or b = 13.08 d, P prec = 162.375 ± 0.011 d, velocity of the ejec¬ 
tions v = 0.2647 ± 0.0008 c, inclination of the jet axis to the 
line of sight i = 78.05° ± 0.5°, and opening angle of the pre¬ 
cession cone 6 = 20.93° ± 0.08°. The ejections and their pre- 
cessi on movement are indeed obs erve d in the radio images: se e 
e.g. Hiellming & Johnston ( 1981 ) and 1 Vermeulen et al Jteb- 
The source has always been observed ejecting material, mean¬ 
ing that this is not only the first microquasar discovered but 
also the only one with continuous ejections, and the only ob¬ 
ject known to show evidence of ions accelerated to relativistic 
velocities (0.26 c). 


Despite intensive studies, the nature of the two stars in 
this binary system remains uncertain. No emission line has 
been identified as clearly belonging to the donor star of this 
X-ray binary, some of them being suspected to be formed 
close to the compact object. Numerous mass ratios were es- 
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timated leading to a low-mass or high-mass X-ray binary 
with either a neutron star or a black hole candidate. The 
presence of Wolf-Rayet-like lines added to the very lumi¬ 
nous continuum in the visible and near-IR ranges led sev¬ 
eral people to associate the donor star to a Wolf-Rayet 
or Of sta r dMurdin et al. 198(1 Ivan den Heuvel et all fl980: 
iHut & van den Heuvelll98ll . However, the coordinated optical 
and X-ray light curves of SS 433 show that the X-ray eclipse 
corresponds to the minimum visual magnitude dStewart et al.l 
1987 1: thus a thick accretion disc is inferred, which is more lu¬ 
minous in the visible range than the donor star. Only recently, 
iChes et ahl ( 2002a ) detected absorption features in the blue 
spectrum suggesting that the donor star is an A-type evolved 
star. In the more detailed spectroscopic study of lHillwig et al.l 
( l2004j) . the donor star is consistent with an A3-71 type with 
M = 10.9 + 3.1 M 0 , and the compact object is a low-mass black 
hole candidate with Mx = 2.9 ± 0.7 M Q . 

The distance of SS 433 has also remained uncertain for a 
long time. From complete observations of the radio emitting 
materials ejected in the relati vistic jets and after taking the kine¬ 
matic model into account, iBlundell & Bowled ( 2004 ) eventu¬ 
ally determined this distance as d - 5.5 ± 0.2 kpc. 

As the origin of the spectrum (lines and continuum) of 
SS433 in the visible and near-infrared ranges appears quite 
complicated for understanding clearly, the idea was to look for 
characteristic lines in the mid-IR range which could help to un¬ 
derstand the nature and the origin of the emissions better and 
possibly to constrain the nature of the donor star further. In 
this article we present observations of SS 433 with the Infrared 
Space Observatory (ISO) in both spectroscopic and photomet¬ 
ric modes in the 2-12/im range at different epochs and at 25 /mi 
& 60 /on. We compare the emission line spectra to the spectra 
of four Wolf-Rayet star of WN type (Sect. 3). We also study 
the continuum emission of SS 433 in Sect. 4 and calculate the 
mass loss and radius of the corresponding free-free emission. 
In Sect. 5 we discuss the constraints on the nature of the mass 
donor star provided by our mid-IR observations, we describe 
the phenomenon imitating a Wolf-Rayet star, and we consider 
the possible large scales behaviour of the strong mass outflow. 


2. Observations and data reduction 

We searched for observations of SS 433, R.A. = 19 h l l m 49 s .57 
Dec. = +04°58'57".8 (J2 000), in the archive s of the Infrared 
Space Observatory 1 (ISO. iKessler et alll 1996h and found sev¬ 
eral of them in 1996 and 1997 as shown in Tabled The main 
observations were achieved with the spectroscopic mode of 
the photometer on board ISO, ISOPHOT iLemke et al.ll 1996 ) 
at four different epochs. The other observations were carried 
out with photometric modes of ISOPHOT and with another in¬ 
strument of ISO: the infrared camera ISOCAM (I Cesarskv et a l. 

s 

We looked in the archives of the Green Bank Interferometer 
(GBI) and at the quick-look results provided by the All-Sky 
Monitor / Rossi X-ray Timing Explorer (ASM/RXTE) team, 
in order to know the state of SS 433 in radio and X-rays, re¬ 

1 http://www.iso.vilspa.esa.es/ida/ 


spectively, at the time of the observations with ISO. The GBI 
2.25 GHz and 8.3 GHz lightcurves show a normal activity in 
the radio, i.e. only small oscillations with no flare during the 
ISO observations. The X-ray level detected by ASM is low (a 
few counts/s), as usual, and has no flaring event. 

In the visible range the spectmm of SS 433 shows strong 
Ha emission lines and looks like the spectrum of a Wolf-Rayet 
star, so we searched in the ISO archives for observations of 
WR-type stars, particularly those of late WN-type which are 
less evolved and so still emit hydrogen lines. We found four ob- 
servations of well-defin ed WN stars (see table|2j with ISOSWS 
rtde Graauw et al Jfl 996l) the short wavelength spectrometer of 
ISO, which has a far better spectral resolution than ISOPHOT. 

2.1. ISOPHOT wide-band filters 

SS433 was observed at 12, 25, and 60 //m with ISOPHOT 
on 11 April 1997 (see Table Q. The measurements at 12 and 
25 /mi were performed in triangular chopped mode; i.e. the 
line-of-sight was switched periodically between SS433 and two 
background positions placed symmetrically at ±120" from the 
source. At 60 /im, also measured in chopped mode, only one 
background position was selected at 180" from SS433. 

The data reduction was performed using the ISOPHOT 
Interactive Analysis Software Package V9.1 (PIA, 
(G abriel e t al. 119971) . After corrections for non-linearities 
of the integration ramps, an 8-point signal pattern was created 
from each observation by overplotting and averaging the basic 
blocks of the observation (the repeated bac kground+source 
cycles); for details see lAbraham et al . l2003t) . The signals of 
the patterns were transformed to a standard reset interval of 
1/4 s, and an orbital dependent dark current was subtracted. 
The signals were corrected for non-linearities of the detectors 
by applying the signal linearization corrections, as available 
in March 2001. From the patterns a [source-background] 
difference signal was then extracted. In most cases the 
measured difference signal underestimates the real signal 
due to short term detector transients: therefore, a correction 
was applied for signal loss. This correction is a function of 
the chopping frequency and the measured difference signal. 
The flux calibration of the 12 and 60 fim measurements was 
performed by comparison with the on-board fine calibration 
source (FCS), which was also measured in chopped mode. 
At 25 /mi the detector’s actual sensitivity could be reliably 
predicted from the orbital position of the observation, and an 
orbital dependent default responsivity was applied. Finally 
the derived flux densities were corrected for the finite size of 
the aperture by using the standard correction values. Colour 
corrections were applied for each measurement. 

The flux uncertainty at 25 /mi was estimated to be better 
than 30%. At 12 /mi the selected aperture of 18" was smaller 
than the standard one in which the ISOPHOT calibration pro¬ 
gramme was executed. At 60 /mi the actual responsivity of the 
detector can deviate significantly from an average default fig¬ 
ure, and its value had to be derived from the accompanying 
FCS measurement, thereby introducing an additional uncer¬ 
tainty factor. Due to these issues, the accuracies of the 12 and 
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Table 1 . Log of the observations of SS 433 with ISO. The T represents the precessional phase from iGoranskii et al . ( 1998h . ch p 
and i if are, respectively, the precessional and orbital phase from lKemp et ali ('1 9861) . The flux densities were given by the PHOT 
data reduction, were directly measured from the CAM images, or were roughly estimated from the ISOPHOT spectra (continuum 
level). See the text for the corresponding uncertainties. 


TDT 

civil date and in Julian Day 

T 

®r 

V 

mode 

A range (FWHM) 
and filters 

aperture 
or PFOV 

flux density 
range (mJy) 

35300404 

3 Nov. 1996 

2450 391.195 

0.41 

0.84 

0.11 

PHT40“ 

2.47M-.87 pm (~ 0.04pm) 

24" x 24" 

380-350 

” 

” ” 

” ” 

” 

” 

” 

PHT40 

5.84-11.62 pm (~ 0.1pm) 

24" x 24" 

330-180 

51200601 

11 Apr. 1997 

2450 550.272 

0.39 

0.81 

0.27 

PHT40 

2.47M-.87 pm (~ 0.04pm) 

24" x 24" 

710-590 

” 

” ” 

” ” 

” 

” 

” 

PHT40 

5.84-11.62 pm (~ 0.1pm) 

24" x 24" 

530-280 

51200602 

” ” 

2450 550.298 

” 

” 

” 

PHT03 6 

PI 12: 11,89 +3.255 pm 

18" 

354 

” 

” ” 

” ” 

” 

” 

” 

PHT03 

P2 25: 23.81 ±4.59pm 

52" 

432 

51200603 

” ” 

2450 550.303 

” 

” 

” 

PHT22 C 

C100 60: 60.8 + 11.95 pm 

43" x 43" rf 

390 

51801005 

17 Apr. 1997 

2450 556.016 

0.43 

0.85 

0.71 

PHT40 

2.47M-.87 pm (~ 0.04pm) 

24" x 24" 

350-450 

” 

” ” 

” ” 

” 

” 

” 

PHT40 

5.84-11.62 pm (~ 0.1pm) 

24" x 24" 

320-200 

52402506 

23 Apr. 1997 

2450 561.939 

0.47 

0.89 

0.16 

PHT40 

2.47-4.87 pm (~ 0.04 pm) 

24" x 24" 

230-330 

” 

” ” 

” ” 

” 

” 

” 

PHT40 

5.84-11.62 pm (~ 0.1pm) 

24" x 24" 

260-180 

70900204 

24 Oct. 1997 

2450746.268 

0.60 

0.02 

0.25 

CAM0L 

LW3: 12-18 pm 

3" x 3" f 

124 

” 

” ” 

” ” 

” 

” 

” 

” 

LW7: 8.5-10.7 pm 

” ” 

148 

” 

” ” 

” ” 

” 

” 

” 

” 

LW9: 14-16 pm 

” ” 

123 

” 

” ” 

” ” 

” 

” 

” 

” 

LW6: 7-8.5 pm 

” ” 

207 

” 

” ” 

” ” 

” 

” 

” 

” 

LW8: 10.7-12 pm 

” ” 

139 

” 

” ” 

” ” 

” 

” 

” 

” 

LW4: 5.5-6.5 pm 

” ” 

226 

71100514 

26 Oct. 1997 

2450748.294 

0.61 

0.03 

0.41 

CAM01 

LW2: 5-8.5 pm 

1.5" x 1.5 " 1 

423 


“ spectrophotometry with ISOPHOT 
b multi-filter photometry with ISOPHOT 
c multi-filter photometry with the far-IR camera of ISOPHOT 
d per pixel, giving a total field of view of 130" x 130" 
e raster imaging with ISOCAM 
f per pixel 

60 pm measurements are expected to be more in the 30-50% 
range. 

2.2. ISOPHOT spectral mode 

Low resolution mid-infrared spectra of SS433 were obtained 
on 3 November 1996, 11, 17, and 23 April 1997, (see TableQ 
with the ISOPHOT-S sub-instrument. The spectra covered the 
2.5^4.9 pm and 5.9-11.7 pm wavelength ranges with spectral 
resolutions of about 0.04pm and 0.1pm and sensitivities of 
26mJy and lOmJy, respectively. The observations were per¬ 
formed in triangular chopped mode with two background posi¬ 
tions located at +90". The orientation of the chopping direction 
on the sky depends on the date of observation and is different 
for each spectrum. The dwelling time per chopper position was 
128 s, and each measurement consisted of 4 OFF1-ON-OFF2- 
ON cycles and lasted 2048 s. 

The ISOPHOT-S data were reduced in three steps:_ 

(1) We used the Phot Interactive Analysis (PIA, iGabriel et al.1 
119971) software V9.1 to filter out cosmic glitches in the raw data 
and to determine signals (in V/s) by performing linear fits to the 
integration ramps. After a second de-glitching step, performed 
on the signals, a dark current value appropriate to the orbital 
position of the individual signal was subtracted. Then we av¬ 
eraged all non-discarded signals to derive a signal per chopper 
step. 


(2) An [ON-OFF] signal averaged for the whole measure¬ 
ment was determined by applying a 1 -dimensional Fast Fourier 
Transformation algorithm; for the application of FFT methods 
for ISOPHOT data reduction see lHaas et alJ (2000)). 

(3) The [ON-OFF] difference signals were calibrated by apply¬ 
ing a signal dependent spectral response function dedicated to 
chopped PHT-S observations iLaureiis et all200ll) . 

In order to verify our data reduction scheme (which is not 
completely standard due to the application of the FFT algo¬ 
rithm) and to estimate the level of the calibration uncertain¬ 
ties, we reduced HD 184400, an ISOPHOT standard star ob¬ 
served in a similar way as SS433. The results are very con¬ 
sistent with the model prediction for the star, and we estimate 
that the systematic uncertainty of our calibration is not greater 
than 10%. However, as PHOT-S was designed for broader spec¬ 
tral features than the lines observed in this study, the measured 
line flux might be underestimated because the ISOPHOT-S pix¬ 
els are separated by small gaps; according to Table 4.5 of the 
ISOPHOT Handbook 2 about 75% of the line flux is measured 
if the line is centred on a pixel’s wavelength; and about 43% if 
the line is centred between two pixels’ wavelengths. 


2 http://www.iso.vilspa.esa.es/manuals/HANDBOOK/pht_hb/node3.html 
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Table 2. Main properties of the WR stars to be compared with SS 433 


Name 

WR type" 

binary 

distance 

A v b 

references 1. 

WR78 = HD 151932 

WN7h WNL 

no 

2.0 kpc 

1.48-1.87 

vdHOl ; C95a 

WR 134 = HD 191765 

WN6b WNE-s 

possible 

1.74 / 2.1 kpc 

1.22-1.99 

vdHOl / M99 ; CS96 

WR 136 = HD 192163 

WN6b(h) WNE-s 

possible 

1.26/ 1.8 kpc 

1.35-2.25 

vdHOl / SH99 ; C95b 

WR 147 

WN8(h) WNL 

B0.5 V at 0.55" 

630 ± 70 pc 

11.6/11.2 

vdHOl / MorOO 


from 
from 
vdHOl = 


Smith etalJll99filt 


van 


der Huchl 1200 lp except for WR 147 


der Huchd l200lli : C95a =iCrowther et alJJl995rl) ; M99 =lMorel et alJll9991) ; CS96 = iCrowther & SmitfJ J 1 9961) ; C95b = 


ICrowther et ~1 J 11995h|) : SH99^ istevens & Howard 119991) ; MorOO = lMorris et ali l2000t) 


2.3. ISOCAM wide-band filters 

SS 433 was observed with the camera of ISO, with 
several large band filters on 24 and 26 October 
1997 (see Table Q. The ISOCAM data were reduced 
with the Cam Inte racti ve Analysis^ software (CIA, 
iDelanav & ISOCAM Interactive Analysis Teamll200db v ersion 

4.0, following the standard processing outlined in lStarck et ali 

(1999). First a dark correction was applied, then a de-glitching 
to remove cosmic ray hits, followed by a transient correction 
to take memory effects into account. Pixels showing remnants 
of these effects were masked, as well as insufficiently lit side 
pixels. The jitter correction was applied for the 1.5" x 1.5" 
resolution image. The flat-field correction used the automatic 
evaluation. Then individual images were combined into the 
final raster map, and finally the pixel values were converted 
into milli-Jansky flux densities. No colour correction was 
applied. We estimated to ~ 20% the total uncertainty of the 
ISOCAM observed flux densities. 

2.4. ISOSWS data on WR stars 

We took the SWS spectra of the four WR stars from the archive; 
i.e. they were processed by the SWS pipeline (officially Off- 
Line Processing = OLP) version 9. In order to compare them 
to the PHOT-S spectra of SS 433, we had to degrade their res¬ 
olution and keep only the wavelength range of the PHOT-S 
modes. The convolution to the PHOT-S spectral resolution was 
performed with IDL routines available on the VILSPA home 
page 3 . 

3. The emission lines: comparison with Wolf-Rayet 
stars 

As we summarized in the introduction, the nature of the visible 
and near-IR spectrum of SS 433 is difficult to interpret, and till 
recently no line could be used to determine the nature of the 
donor star in this X-ray binary. ISO provided unique spectra 
of this object in the mid-IR range. Since SS 433 has for a long 
time been strongly suspected to be a Wolf-Rayet star of WN 
type and its optical spectrum indeed resembles those of WN 
stars, we decided to compare its emission line spectra to that of 
several WN stars found in the ISO archives and rebinned to the 

3 http://www.iso.vilspa.esa.es/instr/SWS/ “SWS vs PHT-S 
Convolution Tool” written by B. Schulz 


I <1> D _ 3 TH_ 

o «-x “o ?~x a> x a) a)J5 ix 

^ cl+ Si Si + ^ + ^+ ++ 

+ ++ + + + ++ TT 

‘2-a57Da5a>®a> CQ -7D 0 — 

mxxxxmxolx Si 



Fig. 1. Observed (i.e. absorption not corrected) spectra of 
SS 433 with ISOPHOT-S on 11 April 1997 and of four Wolf- 
Rayet stars with ISOSWS rebinned to the PHOT-S resolution 
and wavelength range. The A axis has a logarithmic scale, in 
order to better separate the different lines in the low A part of 
the spectra. 

same spectral resolution. The main properties of these four WR 
are shown in Tabled and their observed spectra are displayed 
in Fig. |U along with the spectrum of SS 433 as observed on 
11 April 1997, since it has the strongest emission among the 
available spectra of this source. 

Figure ^ shows the observed spectra i.e. with no correc¬ 
tion for absorption by the interstellar medium. All the emis- 
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Fig. 2. Dereddened (i.e. absorption corrected) spectr a of 
WR 147 and SS 433, using the law of It, utz et all (! 1 9961 ) and 
Ay = 11.2 & Ay = 8, respectively. O d and <£ are the preces- 
sional and orbital phases from lKemp et al.l f 1986 ). as described 
in the text (Sect. IQ . In the spectra of SS 433, some absurd 
points due to instrumental defaults are enhanced by the dered¬ 
dening, as at the upper edge of the 2-5 /mi range, and thus they 
have very large error bars going off of the graphics. 


WR 147 IMorris et alj|200 0). The spectra of SS 433 show the 
same general shape as the one of WR 147 with the same lines 
and comparable relative intensities between the latters, except 
for the 10.5//m line, which is absent in SS 433. Smith & Houckl 
(feooil) studied 8-13 /mi spectra of several WN8 and WN9 
stars, where this [Siv] 10.5//m line can be either very weak 
or not detectable and may be drowned in the noise of the spec¬ 
trum. However, the spectra of SS 433 are not very noisy around 
10.5//m, and the sensitivity of the detector (~ lOmJy) should 
enable us to detect this line if present. The absence of this line 
is surprising, if we are dealing with a strongly ionized, dense 
outflow in SS433, which has a common physical mechanism 
to WN stars. Indeed, one would expect the presence of mid- 
IR fine structure lines, such as [Siv], unless the ionization in 
the relevant region (around a density of 10 5 to 10 6 cuT 3 ) dif¬ 
fers from that of WR stars. Given the similar ionization in the 
He line spectrum, one might suspect that the ionization of the 
plasma is comparable. Thus in the mid-IR, SS 433 (only) looks 
like a late WN star (WNL) of WN8 type, which is relatively 
H-rich (H/He ~ 0.5 - 2. ICrowthetfe003b| ). Note that WN8 stars 
are known to be the most active subclass of Wolf-Rayet stars 
because of their high levels of spectroscopic and photometric 
variability (see iMorris et al.t200oil and references therein). Note 
also that WR 147 is a Wolf-Rayet binary (WN84-B0.5V) sys¬ 
tem with strong colliding wind that has been observed in radio 
(t Williams et al.il 997 


The measurement of line widths and position are lim¬ 
ited by the spectral resolution of the spectrocopic mode of 
ISOPHOT i.e. FWHM ~ 0.04 & 0.1 /mi for the 2-5/im and 
6-12/rm ranges, respectively. Thus it is not possible to detect 
any Doppler shift due to the orbital movement. Concerning the 
0.26 c relativistic motion in the jets, there is no evidence of 
such a Doppler shift in our spectra when comparing the dif¬ 
ferent precesional phases observed, so there is no indication of 
infrared emission from the jets. Because of the limited spec¬ 
tral resolution and the blended lines, we could neither study 
the properties of the lines (equivalent width, shape, etc...) nor 
measure the elemental abundances. 


sion lines were identified using IMorris et Hi] l2000i) and only 
the lines present in at least three spectra of SS 433 and at 2.5 cr 
above the continuum were taken into account. No absorption 
line was found. The emission lines are mainly H 1 lines blended 
with He 1 or He 11 lines. There is also a weak He 1 line at 4.30 fin 1 
and two weak He 11 lines at 3.09 fim and 4.765 /on, but no metal¬ 
lic line is detected. It is clear in Fig. ^ that the spectrum of 
SS 433 is closest to the one of WR 147, a WN8(h) star, although 
the isolated Hei (4.30/tm) and Hen (3.09/mi) lines are more 
marked in this star than in SS433; conversely, the other WR 
stars differ by a strong He 11 3.09//m line. 

The similarities between the spectrum of WR 147 and the 
four spectra of SS 433 also appear clearly once the spectra are 
dereddened, as in Fig. [2] where we corrected the interstellar ex¬ 
tinction using the law front Lutz et a l. (1 9961) (see the next sec¬ 
tion concerning the choice of this law) and a visual absorption 
of Ay = 8 for SS433 ( Mareon 1984) and of Ay = 11.2 for 


4. The continuum 

4.1. Dereddening 

In order to study the continuum emission of SS 433, we have 
to deredden the observed values, i.e. to correct them from the 
absorption along the line of sight toward the source. What we 
call dereddened flux Fd er is linked to the observed flux F 0 b s by 
the formula: Fder = Fobs x 10 +0 4A - 1 , where A.i is the extinction 
at wavelength A, calculated as Ay times the extinction coeffi¬ 
cient at A given by an extinction law (see Fig. 0. As SS 433, 
with l = 39°41'38.8” and b = -02°14'40.6”, is located in the 
galactic plane at a distance of 5.5 kpc, this absorption is not 
negligible. The problem is that the visible absorption Ay is not 
constrained very well for SS 433. 

By fitting the optical continuum slope from UV to near-IR 
with the assumption of various thermal models (blackbody un¬ 
derlying spectral types or/and free-free emission) or using in¬ 
frared spectrophotometry, many authors found Ay ~ 8 (see ref- 
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Absorption laws 



Fig. 3. Compa r ison o f the absorption law bv lDrain IdU with 
the lLutz et alJ dl996l) law in the mid-IR. The coefficient of ex¬ 
tinction is defined by A^/Ay, where Ay is the extinction in the 
optical V (0.55 pm) band and A,\ the extinction at wavelength A. 
The wavelength ranges of the ISOPHOT and ISOCAM filters 
used in this article are overplotted. 


erences in Mur din et alJll98d and lMargonlfl984l) . Evaluations 
using recent X-ray date are not restricting enough, as demon¬ 
strated in [Marshall et al : I2002S) on their model fitting of the 
Chandra grating X-ray spectra of SS433. If they use a sim¬ 
ple power-law model (a poor lit to the data), then they get 
Ah = 9.5 x 10 21 c m~ 2 , which implies Ay = 5.3 using the 
Pred ehl & Schmittl d 1 995l) formula: Ah (cm 2 )/A v = 1.79 x 
10 21 . When they use the X-ray line ratios and simultaneously fit 
Ah, they get 2.2 xlO 22 cm -2 , which implies Ay = 12.29; hence, 
the range of values one will find in the literature. Also from 
Chandra observations. [Namiki et al* (2003) just fitted a power- 
law + gaussian lines and found Ah = (1.0 — 1.3) x 10 22 cm~ 2 , 
which gives Ay = 6.1 - 7.3 mag. In the following we will only 
show spectra dereddened with Ay = 8. Note that using lower 
values such as Ay = 5.3 doesn’t change our results significantly 
(see Sect. 14.41 . 

The extinction of the interstellar medium is highly wave¬ 
length dependent, even in the mid-IR range, but there are dif¬ 
ferent absorption laws available depending on the galactic line 
of sight. The law described by IPrainej (1 9891 ) is commonly 
used, except for regions toward the Galactic Centre for which 
[Lutz et al.l d 19961) showed that a flatter law (between 2.5 and 
8.5 pm) has to be used. Figure l^shows the difference between 
these two laws and the wavelength range of our ISO observa¬ 
tions. The Lutz law is also used in some particular cases like 
some dense H 11 regions or some objects with a carbon-rich en¬ 
vironment. As an example. iMoneti et al.l ( 20oH ) used the Lutz 
law to deredden their 2.5-8 pm ISO data on the Quintuplet 
Cluster, althou gh this young cluster is in the vicinity of the 
Galactic Centre. Ilndebetouw et al.l ( 2005 ) presented very recent 
extinction measurements in the 1.25-8. 0pm range along the 
l = 42° and Z = 284° line of sight in the galactic plane, which 
are_similar and consistent with measurements of ll ,ut.z et al .1 
toward the Galactic Centre. This shows that the Lutz law 


could almost be universal or, at least, more common than previ¬ 
ously thought. iPwe kl J 2 OO 4 I) attributes this law to the presence 
of relatively short metallic needles in the interstellar medium 
that could form in the supernova ejecta and in quiescent O- 
rich stellar outflow. SS 433 may have undergone a supernova 
explosion to form its co mpac t objec t, as it is surrounded by 
the W50 radio nebula (tPubner et al.llil998 li). which is possibly 
a supernova remnant. Thus the Lutz law is likely to be valid 
in the case of SS 433, so we used it in this study. Moreover 
when we dereddened the spectra with the Praine law, the con¬ 
tinuum showed a complicated shape with a narrow “bump” in 
the 7-10pm range, which we were unable to fit correctly, even 
when using a lower Ay value. 


4.2. Variability 


Figure @] shows all the ISO observations of SS 433 once dered¬ 
dened with Ay = 8 and using the law of lLutz etalJdl996l) . The 
continuum level varies from one observation to another, and we 
can try to explain these variations using the orbital and preces- 
sional phases. 


In Table [T] cD p and tp are, respectively, the precessional 


and orbital phases according to (Kemp et aH dl986l) . As usu¬ 
ally ip — 0 corresponds to the primary minimum in the visi¬ 
ble light curve, but here the mid-eclipse is the donor inferior 
conjunction, as the region of the com pact object is more lumi¬ 
nous than the donor star. jKemp et al. 1.1984 chose the preces¬ 
sion phase (D p = 0 at the time of light minimum in terms of 
the 162 d Lourier component, corresponding to the time when 
we tend to look “und er” t he disc, bu t many authors, such as 
iGoranskii et al . ( Il998f ) and [Gies et al .( 2002b ) defined the pre¬ 
cessional phase with T* = 0 (corresponding to <D p = 0.5) when 
the emission lines attain their extremum radial velocities, thus 
the time when the jets are closest to our line of sight and the 
disc presents a maximum cross section to us (see Lig. [5}. The 
values obtained when calculating the orbital phase tp using the 
ephemeris from lKemp et akldl 986l) are quite si milar (±0.01) to 
the on es obtained following the ephemeris of IGoranskii et-alj 
( 1998 ). There is a difference of +0.09 in the results obtained 
for the precessional phases using both ephemeris, probably due 
to the different precess ion periods used by the autho rs: 162.5 d 
forjKempetal. jjl98f ) , instea d of 162.15 d for lGoranskii et all 
dl998l) and lGies et alJd2002bl) . 


We reported the values of the orbital and precessional 
phases of SS 433 with Ligs. 3 and 4 of lKemp et ali ( 19861) . re¬ 
spectively showing the mean V-band light curve on the orbital 
period and 162.5 d light curves centred on four orbital phases 
( 1 p = 0.0,0.25,0.50,0.75). The decreasing level with time of 
the continuum in April 1997 (see Lig. is easily explained 
by the decrease of the corresponding levels in the mean orbital 
light curve and in the precessional light curve. The continuum 
level of the spectrum on 3 November 1996 is remarkably sim¬ 
ilar to the one on 17 April 1997 while the mean orbital light 
curve predicts a slightly lower level. But this could be equi¬ 
librated by its precessional phase which is a bit farther to the 
minimum flux than in April 17. 
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SS433 in April 1997 ; ISOPHOT-S ; Av=8 Lutz law SS433 ; ISOPHOT-S and ISOCAM ; Av = 8 Lutz law 




Fig. 4. Fit of the continuum emission of the dereddened spectra of SS 433 with power laws proportional to T 01 (short dashes) 
or d 0 - 3 (dash dot, left side) or /T 2 (dash dot, right side) or/and /T 0 6 (long dashes) and a black body emission with T = 500 K, 
R = 8000 R q at a distance of 5.5 kpc (dash triple dot). In the left figure the solid line fits the emission of 11 April 1997 (black 
body + power law oc A °- 6 ), and in the right figure it fits the emission of 24 October 1997 for A > 9 jim (black body + power law 
oc A- 2 ). 


However, we cannot clearly explain why the 24 
October 1997 level is so low while its orbital phase <p = 0.25 
should correspond to the maximum orbital flux, unless con¬ 
sidering that the effect of the precession phase is dominant, 
since with <F P = 0.02 + 0.05 this date is very close to the mini¬ 
mum flux when the disc is edge on. Then the much higher flux 
two days later, on 26 October 1997 should also be attributed to 
the slight higher precession phase (<l> p = 0.03 ± 0.05), which 
could lead to a significantly higher flux because of geometri¬ 
cal effects. Note also that these observations in October 1997 
were performed with a different instrument, ISOCAM instead 
of ISOPHOT, so with different sensitivity and spatial resolu¬ 
tion. 



4.3. Data fitting 

The continuum of SS 433 in April 1997 can be fitted very well 
by one or two power laws and a blackbody emission as shown 
in Fig.@](left). For 11 April 1997 the continuum corresponds to 
a power law from 2.5 to 12/rm: F v = C A~ a Jy, where C = 2.3 
and-ar = -0.6+0.05. Then a blackbody model with T = 150K 
and R = 8000 R 0 is added to this power law emission to fit the 
far-IR emission, but with only three points at 12,25, and 60/rm, 
the constraints are not strong. The blackbody model was the 
simplest model to correspond to a possible thermal emission 
in this wavelength range. The two other April 1997 spectra are 
fitted by two power laws with C = 0.6, -a = +0.1 ± 0.15 on 
April 17 and C = 0.35, -a = +0.3^q 2 5 on April 23 for the 


Fig. 5. Convention of the precessional phase O p according 
to |Kem£jjtai]_|i98i^ and the usual precessional phase V F 


(e.g. iGoranskii et al 


shown on the kinematic model of 
SS 433: P or b = 13.08d, P prec = 162.375 ± 0.011 d, velocity of 
the ejections v = 0.2647 ± 0.0008 c, inclination of the jet axis 
to the line of sight i = 78.05° ± 0.5° and opening angle of the 
precession cone 6 = 20.93° ± 0.08° JEikenberrv et alJ200lh . 


2—4.5 jim range, and C = 1.6, -a = -0.6^^ and C = 1.3, 
—a = -0.6^0 j 5 , respectively, for the 4.5-12 /jm part. 

These power laws can be interpreted as free-free emission: 
optically thin for the positive slopes between 2 and 4.5 /vm in 
April 17 and 23, and the negative slope (A~ 0 ' 6 ) corresponding 
to the intermediate regime between optically thin and optically 
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thick free-free emission. This -0.6 slope is the exact theoret¬ 
ical spectral index for the free-free emission from an ionized 
homogeneous wind flowing out from a star with a sp herical 
expansion and at a constant velocity ( Panagia & Fellil ll975 : 
Wri ght & Barlow! 19751) . This is chara cteristic of hot star s with 
stro ng winds, such as O an d WR stars JCohen et al.ll 975|) . Note 
that (Schmid-Burgkl ill 982l) demonstrated that the overall spec¬ 
tral emitting behaviour is the same for non-spherical outflows 
with more complex geometries, as long as they stay thick, i.e. 
when no ratio of structural length scales in the source exceeds 
about 10, which is very likely the case in SS 433. 

Figure 0] (right) shows our fit of the continuum of the 3 
November 1996 spectrum with the same functions as on 17 
April 1997, as they have the same flux level: F v = C A~ a Jy 
with C = 0.6, -a = +0.1 +0.1 and C — 1.6, -a = -0.6 + 0.15. 
These functions fit the continuum well except for A > 10 //m, 
where the continuum decreases more rapidly as if it turns to 
an optically thicker emission. The same behaviour is observed 
with ISOCAM on 24 October 1997, although the wide band 
filters have large error bars. This spectral energy distribution 
follows the Fy — 1.1 A 06 Jy law between 5 and 10/im, then the 
flux decreases indicating a clear cut off toward a larger slope. 
The continuum can then be fitted by an optically thick free-free 
emission (F v = 24 A~ 2 Jy) plus the same black body emission 
( T = 150 K and R = 8000 R Q ) as in the April 1997 fitting. 

Thus over five epochs of observations, the 2-12/rm spec¬ 
trum of SS 433 corresponds to optically thin and thick free- 
free emission, so it is consistent with a standard O or WR 
wind. The far-IR emission is probably thermal emission from 
dust at T = 150 K surrounding the system at a large distance 
(R > 8000 R q ). 


4.4. Mass loss evaluation 


From our fit of the free-free emission of SS 433 and assum¬ 
ing a thick geometry for the wind, we can calculate the corre¬ 
sponding mass loss M (in M 0 yr 1 ) using the lWright & Barlowl 
(1975) equation (8), which gives the flux density F v of the free- 
free emission of a radiatively driven wind: 


F v = 


2.32 x 10 

^kpc 


10 


M 


4/3 


(+GHz y e gZ 2 ) 2li 

RVo, 


mJy 


where 7) kpc is the distance in kpc, ji mean atomic weight per nu¬ 
cleon, v m the terminal velocity of the wind measured in km s _1 , 
vghz the frequency of observation in GHz, y c the number of 
free electrons per nucleon, g the Gaunt factor and Z the mean 
ionic charge. Then the mass loss is 

M= 1,69.10 ~ & F 3 v /4 Dl , p 2 cf iv 00 (v GHz j e gZ 2 r l/2 M 0 yr -1 . 

We took D - 5.5 kpc and we measured the flux densities at 
a wavelength of 5 yum (v = 60000 GHz) in the maximum and 
minimum levels of April 1997: F v = 876 mJy & 495 mJy. The 
other parameters of this equation are not very well known, but 
we can take either the usual or average values of the WN winds, 
assuming that the Gaunt factor isj? =_1 and that the terminal 
velocity is Voo = 1000 km s' 1 dCrowtherl2003ai) . and using a 
wind composition typical of late WN stars /j. = 2, y e = 1, and 


Z = 1 ILeitherer et al ii 1997 ). Note that the adopted velocity 
is also compatible with the observations of th e H„ s tationary 
lines in the optical spectrum; Ivan den Heuvell dl 98 lit reported 
Dopple r widths of ~ 1 - 2 x 10 3 km s _1 , and recent! v lGies et all 
( 2002b ) measured FWHM of ~ 15A (-hSOkms- 1 ). We find 
a maximum mass loss rate of M = 2.9 x 10~ 4 M 0 yr~' and a 
minimum of M = 1.9 x 10~ 4 M 0 yr -1 . 

We explored the possible range of the parameters to see 
how they influence the resulting mass loss. The velocity of 
the wind flowing out from a WR star generally lies between 
500kms -1 and 1500kms , which either reduces the result 
by a factor of 2 or multiplies it by 1.5. In the mid-IR, the 
Gaunt factor is not easy to calculate, as the param eters of 
the plasma are unknown. According to Rvhicki & Lightmanl 
( 1979 ). g ranges from 1 to 5, which would result in a mass loss 
multiplied by 0.45. Another unknown is the chemical compo¬ 
sition of the wind, but the typical WN co mposit ion seems rea¬ 
sonable since the study of ILeitherer et al . (19959 shows values 
of /j between 1.5 and 2.7 for WN7/WN8 stars (which implies 
M x 0.75 to 1.35); and from the mid-IR spectra, the wind is 
dominated by H + and He + ions so y e = Z = 1. We checked the 
results obtained by dereddening the spectra with a lower visi¬ 
ble absorption Ay = 5.3 and found that this does not change the 
mass loss a lot, and it is then multiplied by a factor of 0.85. We 
also calculated the mass loss using the lowest value of distance 
D - 3.2 kpc, inferred from H i observations of the W50 nebula 
around SS 433 JDubner et a lJl998t) . which divides the resulting 
mass loss values by a factor of 2.2. Thus the greatest uncertain¬ 
ties concerning the mass loss come from the wind velocity and 
the Gaunt factor, so that the range 2 - 3 x 10 4 M 0 yr -1 is valid 
within a factor of 3. 

This result agree s well with the past estimates. 
Ivan den Heuvell ( 1981 9 found a wind mass loss rate of 
10 s - 10 4 M 0 yr 1 outflowing from an O or WR type star 
with a velocity of 1000 - 2000km s~'. IShklovskiil ( 1981 ) 
calculated a mass loss rate of ~ 10~ 4 M 0 yr~' from a 10yum 
flux density and by assuming an optically thick plasma with 
a dispersion velocity of v ~ 2000km s~ 3 . More recently, 
iKing et all ( 2000 ) considered possible evolutionary models for 
SS 433 and predicted mass loss rates between 7 x 10 6 M 0 yr~' 
and 4 x 10" 4 M 0 yr _1 . 

We can also compare this result with our previous study of 
CygnusX-3 iKoch-Miramond et al.il2002l) . which is generally 
accepted as an X-ray binary containing a_Wolf-Rayet star. The 
only other such case is IC 10X-1 (Iciark & Crowtheil2004l) . In 
this similar study, the continuum of the mid-IR spectrum of 
Cygnus X-3 is also due to free-free emission and corresponds 
to a mass loss rate of 1.2 x 10 4 M e yr _1 when using the same 
typical WN composition and Gaunt factor and with a velocity 
of Vo o = 1500 km s- 1 . The similarity of both results is striking, 
whereas the observed flux is fairly different, since Cygnus X-3 
is at a distance of 10 kpc with a highly absorbed emission 
(Ay = 20). This shows that the spectrum of SS 433 resembles 
the one of a Wolf-Rayet X-ray binary very well and suggests a 
common physical mechanism at work in both systems. Indeed 
a disc-like wind, though significantly flattened in the plane of 
the binary orbit, has been invoked by iFender et al J (1 999 ) to 
explain both the variability of the He i emission lines in this 
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object and an asymmetric emitting region. This binary system 
has a much smaller size (~ 5 R e separation, 4.8 h period) than 
SS433(13d period), thus the IR emitting region is significantly 
larger than the binary separation, consequently the outflowing 
wind is external to the binary orbit. 

Finally, the mass loss rate of SS433 has to be corrected 
from the effect due to the very probable inhomogeneity of the 
wind. The clump mass loss rates of WN stars are a factor 
of primes lower than previous homogeneous rates dCrowtheti 
l2003a !). so the mass loss rate of SS433 is within the range 
6 - 10 x ICC 5 M q yr _I . This result can be compared to the 
mass lossjate 1.5 — 3.7 x 10 _5 M® yr~' of WR 147 measured by 
iMorris et alJ d2000l) and to the typical range of WN clumped 
mass-loss rat es: 10~ 5 - 5 - 10~ 45 = 3.16 x ICC 6 - 3.16 x 
1O~ 5 M 0 yr _1 JCrowthetl2003bi) . Thus the mass loss rate found 
for SS 433 is the same order of magnitude as a strong WN wind. 

4.5. Size of the emitting region 

Although we don’t know the geometry of this wind, we can cal¬ 
culate the radius of the equivalent spherical free-free emitting 
region thanks to Eq. (11) of Wright & Barlow f 1975 1: 

/ M \ 2/3 

R(v) = 2.8 x 10 28 y l J 3 g 1/3 Z 2 ' 3 T~ 1/2 — cm 

X/dUooV) 

where T e is the electron temperature of the gas (in Kelvin), 
which is as sumed to be consta nt. Combining this equation with 
Eq. (8) of IWright & Barlowl f 1 97 (see Sect. 14.4t gives the 
characteristic radius of the emitting region at a frequency v: 

R(v) = 1.84 x 10 17 T e - 1/2 D kpc v G ^ z < /2 Jy cm. 

The temperature of the plasma responsible for the observed 
free-free emission is unknown, so we take the typical value 
of T e = 10000K for WR stars according to iLeitherer et al.l 
d 19951) . At 5 pm = 6 x 10 4 GHz with the minimum (495 mJy) 
and maximum (876 mJy) flux densities of April 1997, we find 
Rsfim = 3.8 - 5.0 x 10 12 cm = 50 - 67 R 0 . 

The value of this characteristic radius is very dependent on 
the frequency of observation, but also on the plasma tempera¬ 
ture and on the distance. If the plasma temperature varies by a 
factor of 2 or 3, this radius will be multiplied by a factor of 0.7 
or 0.6, respectively. The lower distance 3.2 kpc will multiply 
this result by a factor of 0.6, and the lower visible extinction 
will only multiply it by a factor of 0.9. Moreover, this is the 
equivalent radius referring to a spherical case, so these values 
can only be taken into account as orders of magnitude for the 
typical size of the free-free emitting region in SS 433. 

This result_can be compared to the previous study of 
IShklovskiil dl98ll) . who used a measurement at 10 pm to find 
that the linear radius of the emitting region (considered as op¬ 
tically thick) is R = 5 x 10 12 cm (72R Q ), which is compatible 
with our result. He corrected this with the effects of inclina¬ 
tion of the system to the line of sight, considering the emit¬ 
ting region as an oblate spheroid whose principal radius is thus 
10 13 cm - 144 R a , thus greater than the binary system and 
forming a common envelope around it. These kinds of mod¬ 
els will be discussed later. 


5. Discussion 

5.1. Nature of the donor star and mid-IR constraints 

In the previous section we showed that both the emission lines 
and the continuum of the mid-IR spectrum of SS 433 are com¬ 
patible with a Wolf-Rayet type for the donor star and prefer¬ 
ably with a l ate WN star (WN8). However , t he recent works 
of lciies et alJ ( 2002 a.) and lHillwig et a III 2004) suggest that the 
donor star is an A3-71 supergiant star with 10.9 + 3.1 M®. This 
result seems to be confirmed by the preliminary study of blue 
spectra of SS 433 by ICharles et al.l ( 2004 ). who find absorp¬ 
tion features typically seen in late B-A spectra. On the other 
hand. lLonez et alJ 120041) use a different method, modeling the 
Chandra X-ray spectrum during eclipse to constrain the radius 
of the mass donor star to be 9.1 ± 1.0R®. They note that for a 
main sequence star, this corresponds to a companion mass of 
29 + 7 M 0 , so an 06-08 V type star with T = 41000-35800K, 
according to iLangj (1 992 ). But such a radius could also corre¬ 
spond to a B4-5 giant star. F or a B5III (R = 8 R Q ), we find 
M = 7 M q and T = 15000 K dLangt 19921) . To sum up, follow¬ 
ing all these results the donor star is not a Wolf-Rayet star, so 
there must be a phenomenon imitating such a star in SS 433. 

Before we discuss this point, let us see which constraints 
on the nature of the donor star we can provide from the mid-IR 
spectra of SS 433. In Fig. ©we overplotted the emissions of the 
possible donor stars above cited to the dereddened ISO spectra. 
We modeled these stars with black body emissions at a distance 
of 5.5 kpc. 

Concerning the AI hypothesis, we chose to fit an A7 su¬ 
pergiant which is less luminous (T — 8150 K from lLan ell 9921) 
than an A3 1 (T = 87 70K), and we took a stellar radius of 
R — 60Rq dLangll992l) . This emission is not compatible with 
the mid-IR emission of SS 433, since in the 2-5 fim range the 
A71 flux density is of the order of the continuum of SS 433 on 
23 April 1997 (Fig. ©left). But this continuum corresponds to 
an optically thin free-free emission, which should then not be 
detected because of the dominant black body emission at this 
orbital phase (tp = 0.16) where the donor star is in the front 
part of the orbit. In order to get a black body emission negligi¬ 
ble compared to this optically thin emission, we had to set the 
radius of the A71 star to R = 30 R. or less. Then this constraint 
is compatible with the result of lHillwig et al.l (I2 0041) . who find 
the radius of the Roche lobe volume for the mass donor star to 
be Ri — 28 ± 2R 0 . Note that this corresponds to a star with 
an intermediate size between type I (supergiant) and type II 
(bright giant) stars (lyenij,l99.5l) . 

Concerning a possible donor star with R = 9.1 ± 1.0R®, 
we used the standard values for the radius R and temperature 
T of a B5III and 06-08 V stars, and Fig. ©(centre and right) 
shows that the corresponding black body emissions are negli¬ 
gible compared to the flux density of the continuum obtained 
with ISOPHOT. Thus both a giant B5 star and a main sequence 
06-08 star are compatible with the spectrum of SS 433. 
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Fig. 6. Superimposition of the ISOPHOT mid-IR emissions of SS433 with the flux from the possible donor stars modeled with 
black body emissions at 5.5 kpc: on the left side an A-type supergiant with T = 8150K and R - 60 (dashes) or 30 R Q (dash triple 
dot), at the centre a B5 giant with T = 15000K and R = 8R Q (dash dot) and on the right side a main sequence 06 or 09 star 
with T = 41000K and R = 10 R 0 (dashes) or T - 41000K and R = 8.5 R 0 (dash triple dot), respectively. 


5.2. A phenomenon imitating a Wolf-Rayet star 

Now, as the mass donor star of SS433 is not a Wolf-Rayet 
star, then there must be a phenomenon imitating the emis¬ 
sion of s_uch a star in this binary system. A long time ago, 
Ivan den Heuveli (1 981 1 suggested that SS433 was a system 
containing an early-type star overflowing its Roche-Lobe and 
a large luminous accretion disc, both embedded in a dense, 
spherically symmetric, outflowing atmosphere resembling that 
of a Wolf-Rayet star. We propose a very similar model where 
this Wolf-Rayet phenomenon comes from the material sur¬ 
rounding the compact object and where the donor star is only 
constrained to be smaller than a supergiant. This material does 
not form a classical thin accretion disc but instead a thick torus 
or a thick envelope like a stellar atmosphere, but not spherical, 
which is ionized by the UV/X-rays emitted in the close vicinity 
of the compact object and expelled by radiation pressure, thus 
imitating the wind of a Wolf-Rayet star. 

There may be a “classical” thin accretion disc, but only very 
close to the compact object and not detectable because of the 
surrounding material. This would explain that there is no evi¬ 
dence for an accretion disc_in the Chandra X-ray spectrum of 
SS433 according to [Marshall et ID §002) and that SS 433 is 
not as strong an X-ray source as the other microquasars or X- 
ray transients. 

Note that lHillwig et al.l d2004h find a binary separation of 
a = 0.26 ± 0.02 AU = 56 ± 4 R Q , which is of the order of the 
typical radius of the free-free emitting region at 5 pm R ^ m = 
50 - 67 R b , which we calculated in Sect. 14.51 So the wind may 
envelop a part of the donor star or the whole binary system, 

envelope”, as already 

Several authors have reported the possibility of a pe¬ 
culiar envelo pe sur rounding SS433 with a disc-like shape. 
iZwitter et al. ( 1991 ) proposed that the asymmetric shape of the 
visible light curves is due to a disc-like outflow of matter from 
the accretion disc. This optically thick “excretion flow” is not 


possibly for ming a very hot “c 


confined to the plane of the accretion disc but has an opening 
angle of ~ 20° centred close to the orbital plane and extends 
up to the orbit of the companion. This analysis is compatible 
with the observations of iKotani et al. ( 1199 6). who explained 
that some X-ray lines are heavily absorbed only in the receding 
jet because of the existence of absorbing gas in a precessing 
plane perpendicular to the jets. However, they do not constrain 
the size of this extended rim of the accretion disc, and their 
evaluation of the mass escape rate (a few 10~ 8 - ICR 7 M G yr _1 ) 
is far too low compared to our result based on a spherical model 
dWrivht & Barlowlll 975ll which is also valid f or a thick disc¬ 
like outflow according to lSchmid-Burgk ( 1982 ). Moreover, our 
observations suggest a thicker envelope with roughly the size 
of the SS 433 binary system, more torus-like than disc-like. 

5.3. Ejected material at larger scale 

With such a huge mass loss for SS433 (~ 8 x 1(R 5 M 0 yr -1 ), 
one could ask what happens to all this ejected material. The 
relativistic jets of SS 433 are well known, and the source is 
continuously ejecting material in this way, but (Marshall et alJ 
(2002) and previous studies estimated the mass loss of these 
jets to ~ 10 7 M g yr 1 . Thus the material ejected in the jets is 
negligible compared to the wind that we observe in the infrared. 

A part of this wind material might form dust, which emis¬ 
sion in the far infrared would explain the excess emission that 
we observe for A > 15 pm in the 11 April 1997 spectrum of 
SS433(Fig.g. We modeled this emission with a black body of 
temperature T = 150 K and radius R = 8000 R e — 5.6xl0 14 cm 
= 37.2 AU, which gives a rough estimate of the temperature 
and distance of the dust region. No WN star is known to form 
dust, so this would be another surprising difference between 
SS 433 and WN stars. This dust could form when the ejected 
material cools down, far enough from the binary system so that 
it is not destroyed by the hot radiation of the WR-like object. 
SS 433 would then resemble some planetary nebulae where the 
central star can be very hot (10 4 -10 5 K), particularly the newly 
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discovered [WN7] star associated with the ring nebula PCG11 
JCohen et alJ 12005). The visible s pectrum of SS 433 clearly 
shows C lines: C iii/N m 44644 Aj|Margon[ 1984J iGies et all 
l2002ah and Cn 447231,7236 A JCiies et alJl2()02hl) . The abun¬ 


dance of C may not be high in the wind of SS 433, but the mass 
loss is so huge that it brings a lot of material to the surrounding 
medium, so the conditions for dust formation might be gathered 
in the outskirts of SS 433. 

Fabrikal (fl99i argues that the gas flows out through the 


outer Lagrangian point L 2 behind the compact star and sprin¬ 
kles a spiral of matter forming a disc-like envelope with an 
opening angle of about 20°. In the course of precession, an ex¬ 
tended disc is formed with a double-cone shape with an open¬ 
ing angle of~ 60°. Its inner radius is 10 14 cm, and it is not 
fully filled by the gas up to a distance of 10 15 cm. According 
to this model, in 10 4 years the external radius of this envelope 
is 1.5 pc. Depending on the distance of the L 2 point to the com¬ 
pact object (so depending on its mass) this analysis is compat¬ 
ible (or not), with our thick wind envelope. 

This scenario is also - in general - compatible with the 
radio observations of equato rial outflows in SS433. In the 
1.6 GHz image of IParagi et al.1 11 9991) . there are two faint ex¬ 
tended emission regions at a distance of 30 mas (138 AU) to 
70 mas (322 AU) to the NE and SW from the position of 
SS433. This detection is confirmed bv iBlundell et alJ d200ll) 
who show smooth low-surface brightness emissions extending 
to > 40 mas perpendicularly to the j et at 1 .5 GHz and 5 GHz. 
More observations analyzed bv IParagi et alJ (* 2002 ) show that 
these radio features extending along the equatorial plane of 
the binary system move away from the central engine on 
timescales of weeks to months. It seems that the appearance 
of these radio emitting regions is changing with time, blob-like 
or smoother, and is related to the precessional cycle. The na¬ 
ture of the emission is unknown, very likely non-thermal and 
maybe optically thin synchrotron from a mixed population of 
relativistic and thermal electrons. 

This equatorial outflow of SS 433 has also recently been 
observed with the VLBA bv iMioduszewski et al . ( 120041) . The 
diffuse emission of the outflow seems to brighten as it moves 
away from SS 433 at a velocity of ~ 5000 - 10000km s -1 . This 
is much higher than the wind velocity that we adopted in our 
study and even much higher than the upper limit for the velocity 
of a WR wind. Thus, either the wind outflowing from SS 433 
has a very unusual high velocity - which we cannot explain and 
which would increase the mass loss rate by nearly an order of 
magnitude - or the wind material (mainly electrons) is acceler¬ 
ated with an unknown (magneto-hydrodynamic?) mechanism 
somewhere between the binary system and the large distance 
(> 100 AU) of the radio emitting regions. 


6. Conclusions 

The mid-IR spectra of SS433 obtained with ISO in 1996 
and 1997 show H and He emission lines similar to the ones 
observed in the spectrum of the Wolf-Rayet star WR 147, a 
WN8(h)+B0.5V binary system with colliding wind. The spec¬ 
trum of SS 433 is thus compatible with the presence of a WN8 
star. 


The 2-12//m continuum emission of SS433 corresponds 
to free-free emission, optically thin or intermediate between 
optically thin and thick, depending on both the wavelength 
and time of observation. At 25/mi and 60/im, the emission 
may be due to dust at T ~ 150/if surrounding the binary 
system at a large distance (> 8000 R Q ). Assuming that the 
free-free emission is emitted by a geometrically thick homo¬ 
geneous wind, we calculated the corresponding mass loss of 
~ 2 - 3 x 10 4 M 0 yr~ 1 . If the wind is clumped, this result is a 
factor of 3 times lower: ~ 6 - 10 x 10 5 M G yr 1 , thus compat¬ 
ible with a strong wind from a WN star. 

However, considering recent results discarding a WR star 
for the nature of the mass donor star in SS 433, we propose 
that the WR-like wind observed in the IR is outflowing from 
an envelope of material enshrouding the region of the compact 
object and maybe also the donor star. This envelope, heated, 
ionized, and expelled by the X-ray emission of the compact 
object, is thus imitating a Wolf-Rayet star. This wind is 
probably the source that provides material which forms the 
possible dust emitting in the far-IR and, at larger distance 
(> 100 AU), the equatorial outflow observed in radio. 
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